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Abstract 

We study the evolution of the hybrid entangled states in a bipartite (ultra) strongly coupled qubit- 
oscillator system. Using the generalized rotating wave approximation the reduced density matrices of the 
qubit and the oscillator are obtained. The reduced density matrix of the oscillator yields the phase space 
quasi probability distributions such as the diagonal P-representation, the Wigner lU-distribution and the 
Husimi Q-function. In the strong coupling regime the Q-function evolves to uniformly separated macro- 
scopically distinct Gaussian peaks representing ‘kitten’ states at certain specihed times that depend on 
multiple time scales present in the interacting system. For the ultra-strong coupling realm a large number 
of interaction-generated modes arise with a complete randomization of their phases. A stochastic averaging 
of the dynamical quantities sets in while leading to the decoherence of the system. The delocalization in the 
phase space of the oscillator is studied by using the Wehrl entropy. The negativity of the lU-distribution, 
while registering its departure from the classical states, allows us to compare the information-theoretic mea¬ 
sures such as the Wehrl entropy with the Wigner entropy. Other features of nonclassicality such as the 
existence of the squeezed states and appearance of negative values of the Mandel parameter are realized dur¬ 
ing the course of evolution of the bipartite system. In the parametric regime studied here these properties 
do not survive after a time-averaging process. 



I Introduction 


The quantum mechanical model that describes a two-level system (qubit) interacting with an electromagnetic 
field mode (oscillator) in a cavity possesses a rich theoretical structure. Its dynamical features, particularly in 
the regime of a weak qubit-oscillator coupling and also for a small detuning between the qubit and the oscillator 
frequencies, are described by the exactly solvable Jaynes-Cummings model [T] that employs the rotating wave 
approximation. Recently, however, the stronger coupling domain, where the rotating wave approximation 
no longer holds, has been experimentally investigated. Various novel realizations such as a nanomechanical 
resonator capacitively coupled to a Cooper-pair box [2], a quantum semiconductor microcavity undergoing 
excitonic transitions [3], a flux-biased superconducting quantum circuit that uses large nonlinear inductance 
of a Josephson junction to produce ultra-strong coupling with a coplanar waveguide resonator [3] fall in this 
category. The superconducting qubits and circuits, in specific, offer much flexibility in selecting the control 
parameters. Consequently, they are regarded as suitable candidates for the quantum simulators fMH| . Moreover, 
integrated hybrid quantum circuits involving atoms, spins, cavity photons, and superconducting qubits with 
nanomechanical resonators have triggered much interest [7]. 

To study the strongly interacting qubit-oscillator system where the Hamiltonian incorporates terms that do 
not preserve the total excitation number, the authors of Refs. have advanced an adiabatic approximation 

scheme that holds in the parametric realm where the oscillator frequency dominates the qubit frequency. The 
resultant separation of the time scales facilitates decoupling of the full Hamiltonian into sectors related to each 
scale, and allows approximate evaluation [8] of the eigenstates of the system. To augment the parametric domain 
so that it includes the exact resonance as well as the large detuning regime, a generalization of the rotating wave 
approximation has been proposed m- This generalization utilizes the basis states obtained via the adiabatic 
limit, and subsequently, the excitation number conserving argument d la the rotating wave approximation is 
applied |10j to the Hamiltonian in the new basis. The resultant block diagonalized Hamiltonian now furnishes 
[lOj energy eigenvalues which are approximately valid for the (ultra) high coupling strength and also for a wide 
range of the detuning parameter. 

The hybrid entangled states of the bipartite coupled qubit-oscillator system provide instances of the entan¬ 
glement of the microscopic atomic states and, say, the photonic Schrodinger cat states that may be regarded 
as macroscopic for reasonably large values of the coherent state amplitude. These states play crucial roles in 
many areas. For instance, they facilitate the non-destructive measurement m of the photon number in a field 
stored in a cavity. Bell inequality tests involving these qubit-held entangled states have been proposed m- 
The hybrid entanglement involving a discrete and a continuous quantum variable may offer advantages such 
as achieving near-deterministic quantum teleportation m- Moreover, the quantum bus |14] framework uti¬ 
lizes the hybrid entanglement where direct qubit-qubit interactions are eliminated while providing a commonly 
coupled continuous mode that mediates among qubits. Realization of the micro-macro entangled states via a 
controllable interaction of a single-mode microwave cavity held with a superconducting charge qubit, and the 
subsequent creation of the superposition of macroscopically distinguishable held modes by virtue of measuring 
the charge states of the qubit have been proposed [TS]. Much current experimental activity [TB] is focused 
towards generation of such hybrid entanglement. Recently the optical hybrid entanglement has been observed 
m by the quantum superposition of non-Gaussian operations on distinct modes. Following a measurement 
based procedure, the hybrid entanglement between two remote nodes residing in the Hilbert spaces of different 
dimensionality has been established [15] • The superconducting circuits have also been used for controllable and 
deterministic generation of complex superposition of states mm- 

Our objective in the present work is as follows. Within the generalized rotating wave approximation scheme 
we analytically study the evolution of an initial hybrid entangled state in a coupled qubit-oscillator system. 
Our study includes both the strong and the ultra-strong coupling domains. Tracing over the complementary 
degree of freedom we obtain the time-evolution of the reduced density matrices of the qubit and the oscillator, 
respectively. The qubit reduced density matrix provides the von Neumann entropy of the system that measures 
the entanglement and the mixedness of the state. The oscillator reduced density matrix, in turn, yields the phase 
space quasi probability distributions [25] such as the diagonal P-representation, the Wigner W-distribution, and 
the Husimi Q-function. The Wehrl entropy [24] constructed via the Q-function measures the delocalization of 
the oscillator in the phase space. In the strong coupling regime where a quadratic approximation to the effective 
interaction described by the Laguerre polynomials holds, a long-range quasi periodic time dependence is visible 
for the Wehrl entropy. This is similar to its behavior in the Kerr-like nonlinear self-interacting photonic models. 
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where it has been observed [25ll28| that an initial coherent state therein evolves, at rational fractions of the 
time period, to the superposition of certain macroscopic coherent states popularly known as ‘kitten’ states. 
While paralleling this phenomenon the current qubit-oscillator interacting model has, however, an important 
distinction that stems from the existence of interaction-generated multiple time scales. Quantum fluctuations 
at shorter time scale that signihes the energy exchange between the qubit and the oscillator cause, at rational 
fractions of the long-range time period, a doubling of the number of coherent state peaks (‘kittens’) while 
the Wehrl entropy evolves from a local minimum to a local maximum of the short-range undulations. A 
frequency modulation of the short time scale oscillation is also observed. In the ultra-strong coupling limit 
the randomization of the phases of a large number of incommensurate quantum modes enforces a statistical 
equilibrium, and a stabilization of the occupation on the phase space sets in. Another object of study is the 
Wigner entropy [21] that probes the quantumness of the mixed states linked with the negativity of the W- 
distribution. Our comparison between the Wigner entropy and the Wehrl entropy reveals the close connection 
of their interrelation with the negativity of the W-distribution. Quantum fluctuations in the strong coupling 
domain induce formations of almost pure squeezed states. Using the Mandel parameter [30] the nonclassicality 
of the photon statistics is also observed. 


II The Hamiltonian and the reduced density matrices 


In natural units {h = 1) the Hamiltonian |51fTU] of the coupled qubit-oscillator system reads 

H = ujo)a + -|- \(Tz{a) a), 


( 2 . 1 ) 


where the field mode of frequency oj is described by the annihilation and creation operators (o, a^|h = a^a), 
and the two-level atom having a transition frequency A is expressed via the spin variables (da,, cr^). The qubit- 
oscillator coupling strength equals A. The Fock states {h|n) = n|n), n = 0,1,...; a |n) = \n — 1), |n) = 

\/n + 1 |n-|-l)} provide the basis for the oscillator, whereas the eigenstates CTa,|±l) = ± |±1) span the space of the 
qubit. As the Hamiltonian (2.1) is not known to be exactly solvable, numerous approximation schemes adapted 


to various ranges of parameters have been advanced. For a small detuning between the oscillator and the qubit 
frequencies, and also for a weak qubit-oscillator coupling, the dynamical behavior of the interacting system is 
accurately described [1] by the rotating wave approximation. In another approach, the adiabatic approximation 
[siinj scheme is found to be appropriate in the large detuning limit (A <C to) as it utilizes the difference between 
the time scales of the slow-moving qubit and that of the fast-moving oscillator. Combining the virtues of these 
two approximation schemes a new procedure has been proposed m that maintains a wide range of validity in 
the regime of large values of both the coupling strength and the detuning of frequencies (A ~ 0(a;), A < w). 
This generalization of the rotating wave approximation employs [lOj . as a suitable intermediate step, a change 
of basis to that introduced by the adiabatic approximation IHlISj- Mimicking an inherent feature of the ordinary 
rotating wave approximation, the generalization advanced in m retains, in the transformed basis, the ‘energy- 
conserving’ one-particle transition elements in the Hamiltonian matrix which now assumes a direct sum of 2 x 2 
diagonal blocks apart from an uncoupled ground state. 

For the purpose of making our notations and subsequent applications clear, we, following |10j . now briefly 
review the diagonalization of the Hamiltonian in the generalized rotating wave approximation scheme, and 
explicitly list its eigenstates. The adiabatic approximation physically signifies that the rapidly moving oscillator 
adjusts itself quickly to the slow changing qubit observable cr^, which may be considered to reside in one of its 
eigenstates | ± 1). Setting A = 0, the bipartite Hamiltonian ( |2.1[) is readily diagonalized in the tensored basis 
|il,n.±), where the displaced number states read: |n±) = (±A) \n) , D (a) = exp (aa^ — a*a) . 

eigenvalues of both the displaced basis states |±l,n±) are degenerate: = |),a: = 4A^/w^ 

elements for the projections of the displaced number states |8] are listed below: 


a € C. The 
The matrix 


Mm,n = {m-\n+) = 


exp ( — f) ^/nXjrn}. L 
exp ( - I) y/m\ln\ Lm~"'\x) 


(m—n) 


(x), m> n 
in < n, 


( 2 . 2 ) 


where the associated Laguerre polynomial reads Ln\x) = X]r=o (n-l) ff- reality of the matrix 


element (2.2) leads to the identity 


Mr. 


= (-l) 


m+n 


Mr 


(2.3) 
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that will be used later to simplify the density matrix elements. In the adiabatic approximation scheme, the 
degeneracy of the states are lifted by the qubit Hamiltonian that causes the mixing of the equi-displaced bipartite 
states |±l,n±). Diagonalizing the resulting 2x2 blocks associated with the same photon number the energy 
eigenvalues and the corresponding eigenstates are obtained [5]: 


p(=t) _, , 

"(>o) “ 


^ ^ = ^exp(-a;/2), |£;W) = ^(|1,n+) ± |-1,n_)). (2.4) 


The basis obtained in (2.4) is now utilized m for presenting the generalized rotating wave approximation 
scheme. The Hamiltonian (2.1), when expressed in the adiabatic eigenstates (2.4), shows characteristics similar 
to the standard rotating wave approximation. Proceeding parallely, the tridiagonal matrix elements that contain 
the ‘energy conserving’ one particle transitions are retained [lOj . and the remote off-diagonal matrix elements 
that manifest higher order transitions are neglected within the generalized rotating wave approximation scheme. 
The ground state is now decoupled from all other states, whereas the remaining matrix elements are organized 
in a tower of 2 x 2 diagonal blocks m as 


E. 


(-D 


Cn 

A-) 


Cn = Y Mn-l,n = ^ , U > 1. 


V Cn K 

The energies of the singlet ground state, and the infinite tower of the doublets read uni 


’(±) 


(2.5) 




^C-) 
— ^0 

tox 

A 

~ 4' 

~ ~2’ 

( 

1 X 

\ 

= (jj[n 



V 

2 4 

) 4 V 

1 [, 

r A / 

L„_i(ai 




A2 


The corresponding eigenstates are explicitly given by 

|£o) = |4”^) = ^( |1,0+) - 1-1,0_)), |£(g,)) = ± 


( 2 . 6 ) 


(2.7) 


where we abbreviate: = ^Cn + ^n; ^n = 

bipartite basis states (2.7) now reads: 




— 


(±) _ / X,^±e 


2Xr. 


^. The completeness requirement of the 


OO OO 

|fo)(4l + E + 14”^)(4”^l) = E + \Ei-^){Ei^\) = I (2.8) 


n—0 


With the above construction of the energy eigenstates via the generalized rotating wave approximation 
procedure in hand, we now study the time evolution of the entanglement of the bipartite system. The initial 
qubit-oscillator hybrid entangled states read 


l'I^(0))(±) = ^ (|l)(|a) + |-a))±|-l)(|a)-|-a))) , |a) = D(a) |0) . 


(2.9) 


As the Schrodinger cat states ^ |a) ± |—a) are mutually orthogonal, the initial states (2.9) may be regarded 
as examples of the hybrid Bell states. Even though our subsequent analysis may be developed with both the 
initial states (2.9), we, for the purpose of notational simplicity, quote the results for the choice |'I'(0))(_j and 
omit the subscript hereafter. The complete basis states (2.7) immediately implement the mode expansion 


|vi/(0))=Co|£:o) + Eci^^ 


( 2 . 10 ) 
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The coefficients of the vector space expansion (2.101 are listed below: 




1 

71 

1 

71 


exp y- 
exp 
xp |- 


f |a_P \ T / ^ 

-^-, <Pa = -Im(a), 


2 

l«+P 


n—1 


(±) _i_ (t) 

f^n ^ /— I > I 

V'^ICnl 


•7 ( 


(±) _ 

r'n ' ! — I /- I t^n 

'Jn ICnl 


)] 


( 2 . 11 ) 


where the following definitions are used: a± = a± • The basis of approximate eigenstates 


(2.7) of the Hamiltonian (2.1) readily yields the time evolution of the initial state (2.10): 

OO 

|vi/(t)) = Co(t) |fo) + E ^710 


( 2 . 12 ) 


n—1 


where the coefficients read: Co(t) = Cg exp (t) = exp ^ . The time dependent density 

matrix of the bipartite pure state has the tensored form 


p(t) = |vl/(t))(d^(t)|. 


(2.13) 


The reduced density matrices of the individual subsystems are extracted by partial tracing over the com¬ 
plementary subspaces in the full Hilbert space. For instance, the tracing over the oscillator states yields the 
reduced density matrix of the qubit: 


_ rp S ^ 

PQ = Tro p = 1 _ 

where the matrix elements read 

OO 

Q = Re^Co{t)Ai{t)* +'^An+lit)*Bnit)^ , 

n—1 

1 ^ OO 

^ = --|Co(t)pMo.o- 2 E ((^oW*-T„(t) + (-l)"Co(t)A(A)*)A^o.„-i 

n—1 

^ 00 

+ (Co(t)*H„(t) + (-l)”Co(A)H„(t)*)Mo.n)+2 E (A^(i)*A„(i)M. 

n,?7i=l 

-H„(t)*H„(t)M„,„ + (^An{t)*B^it) + {-ir+^An{t)B^it)*)Mn-l,m] , 


(2.14) 


In—l.m—1 


An{t) = p7)c(+)(t)+/i7)ci-)(t), B^it) = ^ (2.15) 


In the above expression we have used the identity (2.3). The expectation value of the qubit spin variable (az) 
measures the statistical average of the population inversion. For the density matrix (2.14) it is given by 

(az) = Tr{azPQ) = 2g. (2-16) 


In general, the pair of eigenvalues of the qubit density matrix (2.14) 

i ±ro, G7 = 77 T 77 

allows us to compute its von Neumann entropy S{pq) = —Tr(pQ logpg) as 

5'c = - ( ;^ + ^^ ) log f E w ) - f J - w ) log ( i - 


(2.17) 


(2.18) 
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The bipartite density matrix (2.13) also similarly produces the reduced density matrix of the oscillator via 


partial tracing on the qubit Hilbert space: po = Trg p. Its explicit construction is listed below: 

OO 

n—1 

OO 


-‘rPn(t)’Am(t) Pn^m-1 + P^_{ 


(-) 


(2.19) 


where the symmetrized projection operators read Pn^ = h (l^+) (^+1 i |n_) (m_|), (n, m = 0,1,...). The 


density matrix (2.19) obeys the normalization condition: Trpo(t) = 1. For later use we quote the initial 
oscillator density matrix: po{0) = 5 (la) (a| + |—a) (—a|)- 

To obtain measures of the entanglement and the mixedness of the bipartite system the von Neumann entropy 
of the reduced density matrix of the qubit may be considered. It is well-known m that if a composite system, 
comprising of two subsystems, resides in a pure state, the entropies of both subsystems are equal. In the present 
example it holds for the oscillator with an infinite dimensional Hilbert space, and the two-level qubit interacting 
with it: Sq = So = S. 

Ill Quasi-probability distributions on the oscillator phase space 

III.A The diagonal Sudarshan-Glauber P-representation 

It is well-known that the coherent state representation is overcomplete [23]. Employing the overcompleteness 
property, it has been established |32l |33| that an arbitrary oscillator density matrix po may be realized via 
diagonal coherent state projectors: 

Po = j P{l 3 ,l 3 *)W){l 3 \d^l 3 . (3.1) 

The normalization of the density matrix po ensures the property 


J P(/3,/3*)d2/? = l, 


(3.2) 


whereas its Hermiticity implies P{j3,j3*) to be real. The P-representation, however, is a general distribution 
of indeterminate sign as it possibly includes highly singular derivatives of the (5-function [52] ■ For the coherent 
state, which corresponds as closely as possible to a classical harmonic oscillator state with a complex amplitude, 
the P-representation is a positive definite (5-function [S^ which is a legitimate measure of the classical probability 
density. For this reason the coherent state is considered to be classical. On the other hand, if the diagonal 
kernel P(/3, /?*) can not be considered a valid probability measure due to its negativity, the oscillator state given 
by the density matrix po exhibits nonclassical features. The diagonal P-representation is known to be unique, 
even though it may be realized in different equivalent forms [34j . Towards inverting the diagonal representation 


(3.1) the Fourier transform is employed 


F(/3, P*) = ^ ^ J (- 7 I Po It) expdTp) exp(/ 37 * - /3*'y) d^ 7 . 


(3.3) 


The invertibility of the relations (3.1 3.3) suggests that both the density matrix po and the diagonal rep¬ 


resentation P(P,/3*) encapsulate equivalent informations. We now explicitly evaluate the pseudo probability 
distribution corresponding to the oscillator density matrix (]2.19 ): 


OO 

P{f3,n = -|Co(t)pp(+^+Re(Co(t)*^(A((i)4:li+S„(i)4+J) 

n—1 

00 - 00 

-f H„(<)M™(t)P^_i) + - ^ (A„(i)M™(t)P, 


(-H) 

n—l,m—1 


n,m=l 


n,m=l 




(3.4) 
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where the arbitrarily singular derivatives of the (^-functions read 


= 


1 

V n\m\ 




d 

dj+, 


±exp(|^_n(- J-)>)(/?_)) , /3. = /^± f . 


5/3_ 


(3.5) 


The diagonal P-representation (3.4) satisfies the normalization condition (3.2), and is obviously real. At t = 


0, the initial hybrid entangled state ( |2.9[ ) gives rise to a positive definite P-representation: P(/3,/?*)|(=o = 
2“^ exp(—(|ap — |/?P)) {a —13) + {a + (3)). At an arbitrary time, however, a large number of interaction- 
dependent modes develop, and produce quantum interferences among themselves. The quantum properties 


in (3.4) are evident due to the appearance of rapidly oscillating singular derivatives of 5-function in the P- 


representation, which at time t > 0 is clearly not a nonnegative definite quantity. The existence of negativity 


of the diagonal P-representation (3.4) points towards the manifestations of specific quantum effects such as 
anti-bunching [3S] and quadrature squeezing [3B]. Other phase space pseudo probability functions such as 
Wigner IT-distribution and the Husimi Q-function may be obtained from the diagonal P-representation after 
its suitable smoothing by the Gaussian kernels. 

III.B The Wigner W-distribution 


The Wigner phase space distribution for a quantum state given by the oscillator density matrix (2.19) is 
described [23] as 


W(/3,r) = ^J Tr(D( 7 )po)exp (/Sy* - ^*7)d'7, 
which is a real-valued function maintaining the normalization property: 

I W(^,/?*)d 2 /? = l. 


(3.6) 


(3.7) 


A characteristic of the Wigner distribution that distinguishes it from a classical probability density, is that its 
integral over a given subregion of the phase space may be negative or greater than one. As an efficient compu¬ 
tational algorithm the W-distribution has been expressed m as a sum of an infinite series of the expectation 
values of the oscillator density matrix in the displaced number basis states: 


W(/?, 13*) = - f](-l)" (/3, k\po 1/3, k ), |a, n) = D(a) \n). 

TT • * 


(3.8) 


fe =0 


On the other hand, the convolution of a smoothing Gaussian function of variance 1/2 on the phase space with 
the diagonal P-representation leads to the W-distribution that is free from the divergences brought forward by 


the derivatives of the 5-function contained in the expression (3.4) of the P-representation: 


W{P,n = l J P(7,7*)exp(-2|/3-7nd27. 


(3.9) 


Employing the explicit evaluation (3.4) of the P-representation and the convolution integral (3.9) we now provide 


a direct derivation of the W-distribution. Towards this we employ the identity 

yexp(-2|/3-7p + l7p)(-^) (" ^) 

= 2"+“/3"/3*™ exp(-2|/3n 2 Po(- n, -m; , 


(3.10) 


k-1 


where the hypergeometric sum reads V?-J'^) = {x + £). The above hyperge- 

k=o ■ i=o 

ometric function with the negative integral numerator may be expressed [38j via the Charlier polynomial: 
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c/c(£;r) = 2^0 Vr > 0. To emphasize the equal footing of the integers {k,tj in (3.101 we, how¬ 

ever, do not use this notation here. The integral representation (3.9) aided by the identity (3.10) now generate 
the IT-distribution: 


OO 

w{p,n = \co{t)f +2 Re(co(t)* h) 

n —1 

OO OO 

+^ B^{trA^mi:LiiP, n ] +E n 

n,m—l 


n,m—l 


(3.11) 


where the weights involving the Gaussian functions are given by 

on+m . , 1 . 

[^+P*+ exp(-2|/3+n 2 Jb( ) 

TTVnlm! t ^ ^ V 4|d+|"^/ 

±/3!!/3r exp(-2|/3_n2Fo(-n,-m;_;-^^)) . (3.12) 

An alternate derivation of the Wigner IT-distribution follows from sum rule ( |3.8[ ). Towards this we list the 
matrix elements 


(n±|a,fc) = (-1) 


k '^±^± 

\/n\k\ 


exp 


( 


a± 


T '^4*0^ 2^0 ^ 


- n, -/c; -T 


|a±| 


(3.13) 


that facilitate the evaluation of the series (3.8). Another necessary tool for the present derivation is the identity 


' , 1\ / In 

H - 2 i^o( - n,-k;-;--) 2 Fo(- k, -m; j = 

k^O 

2"+™ exp(-r) 


(3.14) 


that readily follows [38] from the bilinear generating function of the Charlier polynomials. The series sum (3.8), 


in conjunction with the results (3.13 3.14), now reproduce the IT-distribution obtained in (3.11). This serves 
as a consistency check on our derivations (3.4 3.11) of the phase space distributions in the present model. 


The negativity of the Wigner function is an indication towards the existence of nonclassical properties 
of the states. At t = 0 the W-distribution for our initial hybrid entangled state (2.9) is positive definite: 
W{f3, /3*)\t=o = 7r“^ (exp(—2|a — /3p) -I- exp(— 2|q! -I- /3p)). At a later time, however, interactions produce a 
large number of modes that interfere with each other. Interference fringes gives rise to appearance of phase 
space domains representing negative values of the W-distribution. A fruitful measure of negativity relate to the 
volume of the negative part of the Wigner function on the phase space [55] : 


Sw = I |W(/3)|d2/3-l. 


(3.15) 


We will return to the topic in Sec. IV 


III.C The Husimi Q-function 

The Husimi Q-function |23j is a quasi probability distribution defined as expectation value of the oscillator 
density matrix in an arbitrary coherent state. It assumes nonnegative values on the phase space in contrast 
to the other phase space quasi probabilities. Being easily computable it is extensively employed HOI mi in 
the study of the occupation on the phase space. For our reduced density matrix of the oscillator (2.19) the 
Q-function 

Q(Ar) = -(/3|po|/3) (3.16) 
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maintains the normalization restriction: / (5(/3,/3*)d^/3 = 1 and the bounds: 0 < Q{(3,I3*) < Our construc¬ 
tion of the oscillator density matrix (2.19[) now yields the time-evolution of the Q-function: 


n—1 

oo 


72,771—1 


+ 


. OO 

- ^ (3.17) 

72,772 — 1 


where the weight functions on the phase space read 

1 


= ^^(/3+/3r exp(-|/3+n±/3!!/3r exp(-|/3_n). 

'Kynlm: ^ ^ 


(3.18) 


A convolution integral of the diagonal P-representation with a Gaussian weight of unit variance produce 
the well-behaved nonnegative Husimi Q-function as follows: 


Q(/3, 13*) = ^ J P{l, 7*) exp(-|/3 - 7 |^) d^y. 


(3.19) 


The explicit constructions of the P-representation (3.4) and that of the Q-function (3.17) obey the consistency 
check (3.19). Another convolution property that express the Q-function via the Vk-distribution reads [53]: 


Q{(3,(3*) = ^ J ^(7,7*)exp(-2|/3-7p) 


(3.20) 


Employing the integral representation (3.10) it follows that the phase space quasi probability densities (3.11 


3.17) maintain the integral sum rule (3.20). The convolution relations (3.9 3.20) suggest an interesting feature. 


The diagonal P-representation containing singular derivatives of d-functions is smoothed at two stages by a 
Gaussian weight of variance 1/2. After the first smoothing operation the VE-distribution, which is nonsingular 
but not non-negative, appears. Subsequent smoothing operation produces a nonsingular and nonnegative Husimi 
Q-function. If a single smoothing operation (3.19) is performed on P-representation via the Gaussian weight of 
variance 1 , the Q-function is directly reproduced. We symbolically summarize this as [[P] 1 / 2 ] 1/2 = where 

the notation reads [X]s = ^/ <^( 7 , 7 *) exp ^ d^ 7 . 

One of the utilities of the Q'function is that it provides a convenient evaluation of the expectation values of 
the operators expressed in their antinormal ordered form |23| . The first and second moments of the quadrature 
variable defined as 

Xg = ^ (aexp(—i0) + exp{i6)) (3.21) 

may be expressed [23] via the following phase space integrals over the complex plane: 

{Xg) = Tr (Xgpoit)) = \j (^e-*^ + l3*P^)Q{f3,13*)d^l3, 

(3.22) 


(3.23) 


(X^g) = TiiX^gpoit)) = -J ((^e-*^ + /3*P«)2 - ljQ(/3,r)d2^. 

The Q-function (3.17) now allows explicit determination of the above quadrature moments: 
(Xg) = Re(^Gi{t) exp{—i0)^ — ^/xgcosO, 

{Xg) = i Re^G 2 (t) exp(-i26i) - a/x Fi(t) (1-I-exp(-i26»))^ 

+i(N.(i) 


— X I + V cos" 


(3.24) 

























where the time-dependent coefficients {Gfc(t),Nfe(t),F/j(t)|fc G (1,2)} expressed as sums over Fourier modes 
read: 


Gk{t) = Vk Co{t)*Bk{t)+'y'^ \y'{n)k An{t)*An+k{t) 

n—1 

+ \/(n + l)fe Bn{t)*Bn+k{t)^ , 

oo 

k{t) = |Co(t)p-f^ [n'^ An{trAn{t) + {n+lf BnityBnit)], 

n—1 

oo 

Ffc(t) = k Co{t)*A 2 {t) + '^ ^^/n {n + l)k-iAn{t)*Bn{t) 

n—1 

+ \/ri-\-l{n + 2)k-l Bn{t)*An+ 2 {t)'^ ■ 


N 


(3.25) 


The variance of the quadrature variable is given by Ve = (Xg) — {Xgy, which we will later employ for studying 
emergence of the squeezed states during the time evolution. For later utilization we also quote here the mean 
photon number (h) = Tr{Aa po{t)) and its variance ((Ah)^) = Ti{iA po{t)) — (h)^: 


(h) = Ni(t) + ^- y/xRe [Fi(t) j 


- 1 


((Ah)2) = N 2 (t) + |(Ni(t)+Re(G 2 (t))-^) +VS(Re(Fi(t))-2Re(F2(t))) 
- (Ni(t)-v^Re(Fi(t))) 


(3.26) 


(3.27) 


Another dynamical quantity that is useful in the study of ‘kitten’-like states is the polar phase density of 
the Husimi Q-function m obtained via its radial integration on the phase space: 


poo 

Q{e)= Q(/3,r)l^|d|/3|, /3 = |/3|exp(*0), 

Jo 


(3.28) 


which is a convenient tool for describing the splitting of the Q-function. Towards obtaining a series expansion 
of the phase density Q{0) we define the weight factors depending on the polar angle: 


ny;liO)= i7W(/3,/3*)|/3|d|/3| 


(3.29) 


that admit, via (3.181, explicit evaluation as follows: 
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(3.30) 


where we have used the integration given by 

^ rPexp(-(r ± A)2)dr = 


2 ’ 2 ’ 


-A 


0 - 


(3.31) 


The above compendium of expressions (3.28 3.17 3.30) now provides the evolution of the polar phase density 
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Q{9) for our initial state (2.9) 


QiO) = -|Co(t)pH(+)(0) + Re(Co(t)*5](^„(t)H<,7li(0) 
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oo 
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IV Delocalization on the phase space 

IV.A Wehrl entropy 


(3.32) 
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(b) 
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Figure 1: The time evolution of the Wehrl entropy Sq for the parameters A = 0.8 w, a = 2.5, and various 
values of A: (a) 0.008 w, (b) 0.01 w, (c) 0.012 a;, (d) 0.1 a;. In the cases (a)-(c) a quasi periodic behavior 
is observed, whereas in (d) the randomized phase relationship between a large number of modes leads to a 
stochastic stabilization of Sq. 


The Wehrl entropy defined as 


Sq = - Q(/3,/3*) logQ(/3,r)cl2/? 


(4.1) 


is an information-theoretic measure estimating the delocalization of the system in the oscillator phase space. It is 
considered [42j as a count of an equivalent number of widely separated coherent states necessary for covering the 
existing phase space occupation of the coupled oscillator. Being subject to the restriction originating from the 


Heisenberg uncertainty principle, the Wehrl entropy (4.1) is a positive definite quantity [24]. In the present case 
we employ the definition (4.1) and the evolution (3.17) of the Q-function to numerically study the long-range 


time dependence of Sq for various values of the coupling strength. We note that here and hereafter the time is 
measured in the natural unit: u!~^. We observe the following properties: (i) In the long time limit t > (a;^A)“^ 
the quasi periodicity of the Wehrl entropy is manifest in the coupling strength regime A < 0.05 a;, where the 
Laguerre polynomials Ln\x) are well-approximated by their quadratic components 0{x^). Frequency modes 
0{x^A) and their harmonics produced via the interaction now give rise to the quasi periodicity of Sq, where 
the long range time period maintains the property: Tiong ot A”"^ exp(a;/2) (Figs. (a)-(c)). For a smaller value 
of the qubit frequency A the quasi periodic behavior persists for a comparatively higher coupling strength A. 
This follows from the requirement that for the quasi periodicity to hold, the phase change caused by the higher 
order fluctuations {0(a:"’)]ri > 2} during the time span Tiong is to be negligibly small: rionga:^A ^ 1. The time 
period observed in Figs, ^a)-(c) are noted below: Tiong = 6.9914 x 10® (for A = 0.008a;), Tiong = 2.9568 x 10® 
(for A = O.OlOo;), Tiong = 1.4743 x 10® (for A = 0.012a;), respectively. The near equality of the product 
Tiong A^ exp(—a;/2) in the respective cases (0.02863,0.02956,0.03056) validates our argument that the quantum 
fluctuations 0{x'^) produce the observed long range time period. In the instance of nonlinear Kerr-like medium 
similar behavior in the time evolution of Sq was previously noticed |26] . where its local minima corresponded 
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Figure 2: For the parametric values of A = w, A = 0.3 w, a = 4 (a) refers to the time evolution of Sq. At its 
local minimum (wt = 77) the Wigner fF-distribution (b), the Husimi Q-function (c), and the polar plot (w.r.t. 
0 ) of the phase density Q{0) (d) are given, respectively. The localized single peak of the phase density Q{9) 
occurs at 9 = 355.49°. 


with the formations of finite superposition of coherent states. In the present bipartite interacting rnodel, 
however, the qubit-oscillator interaction superimposes short time span fluctuations of the frequency 0{xA) on 
the long range oscillations that now act as an envelope of the total time evolution of the Wehrl entropy. This 
introduces important distinctions to the present model. We will discuss this in the Subsec. IV.B (ii) In the 


ultra-strong coupling regime A > O.lw all frequency modes {0(x"’A)|n = 0,1,...} and their harmonics arise. 
Random phase differences between a large number of incommensurate modes cause the resultant interference 
to average out, while ensuring an effective stabilization of the occupation of the phase space after an initial 
build up (Fig. (d)). The rapid high frequency 0{oj) fluctuations are of small amplitude: |AS'q|/S'q ^ 1, 
and may be removed by a suitable coarse graining process [43) . In this regime it is observed that for a fixed 
initial state parameter a the time-averaged value of the Wehrl entropy gradually increases with increasing 


coupling strength. Following (3.26) it is evident that higher coupling strength leads to an enhancement of the 
photon expectation value (h) that causes a wider spread of the Q-function resulting in an increment in Sq. (iii) 
Sufficiently localized states are transiently observed at large macroscopic values of the coherent state amplitude 
a even in the ultra-strong coupling domain. These states are realized (Fi g© at the local minimum of the Wehrl 
entropy Sq. Physically, the energy exchange between the qubit and the oscillator degrees of freedom induces 
the revival and collapse of the qubit density matrix elements. A revival of the qubit matrix element indicates 
that the oscillator has, reciprocally, less energy available to it. This acts as a constraint on its delocalization on 
the phase space. Moreover, a large value of the amplitude a signifies comparatively higher magnitude of energy 
residing in the interference pattern (Figj^ (b)) of the quantum oscillations. Qualitatively, when the energy 
associated with the stochastic randomization of the modes is less than the energy of the coherent interference 
pattern, a localization on the phase space takes place and, consequently, a relative decrease in Sq develops. 
The quantum interference pattern (Fig[^ (b)) has the shape of a localized peak with ‘twisted arms’, where 
oscillations develop perpendicular to these arms causing a transport of energy necessary for the localization. 
The phase density diagram (Fig. (d)) shows that the initial state (2.9) consisting of two almost maximally 
mixed Gaussian peaks coalesce at the local minimum of Sq to produce a partially pure transient state of the 
oscillator with its von Neumann entropy given by 5" = 0.5883. (iv) The fast initial rise of Sq in the ultra-strong 
coupling domain A ~ w may be understood as follows. Increased qubit-oscillator coupling leads to the generation 
of all high-frequency quantum fluctuation modes. The phase randomization of the modes of incommensurate 
frequencies results in the rapid initial spreading on the phase space, and consequent fast production of Sq. 
Once the modes are statistically populated, the Wehrl entropy S'q, except for high frequency 0{iS) quantum 
fluctuations of relatively small amplitude, maintains an almost stationary value. The initial production time of 


the Wehrl entropy T,„t. prod, follows from the asymptotic behavior of associated Laguerre polynomials L 
large 1, fixed j, and a: > 0 [44] : 


U) 


(x) at 




. n 2 4 exp(a:/2) 

X) = — -j= --■— I — cos 


nx--\j + l 


+ Oin^ 




(4.2) 
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The asymptotic limit of the energy eigenvalues (2.71 is now readily obtained: 


^ ^ ~ \ ~ ^ ^ ^ cos (2v^ - 0 + 0(n ^), (4.3) 

where the leading interaction-generated part o{n~^'^ on the rhs provides the effective statistical stabilization 
of the occupation in the phase space. Recognizing this, the coupling strength dependence of the typical time 
scale for the generation of Sq is now given by Tent.prod. (nx)~^ ~ 1 ^ Tent.prod. -^/A. The rapid initial 

increase of the Wehrl entropy is described in Fig. For the data presented in Fig. the proportionality 
constant Tent.prod./V^ read 30.99,29.66,33.11 for the coupling strengths A = 0.9 w, 1.1 w, 1.3 w, respectively. 
The discrepancy 6%) in the observed data occurs since the local fluctuations in the time evolution of 
play an important role in determining Tent. prod.- A suitable coarse-graining process [43] to smooth the high 
frequency fluctuations may be adopted for fuller agreement. 

Figure 3: The initial rise of the Wehrl entropy Sq for the pa¬ 
rameters A = 0.5 01,0 = 2.5, and the coupling strength A = 
0.9a;(blue), l.la;(green), 1.3a;(red). The horizontal dotted lines 
(for the respective colors) measure the corresponding time-averaged 
values of Sq that read 4.0420,4.1280,4.1890 in the said order. The 
averaging is done for a time interval ojt = 4000 in each case. The 
vertical dotted lines signify the first crossing of the averaged Sq by 
the corresponding time evolution graph. This is regarded as the 
measure of Tent.prod.- The arbitrariness in this determination can 
be improved by a local smoothing operation. For the respective 
coupling strengths these estimates read 29.40,31.11,37.76 in the 
said order. 



Ljt 


IV.B Kitten states and multiple time scales 


For the Kerr-type nonlinear self-interacting photonic models, the local minima in the time evolution of Sq of 
an initial coherent state are associated [26] with transient formations of the superposition of a finite number of 
coherent states [251128] maintaining a uniform angular separation on the complex plane. These superpositions 
are realized at rational sub multiples of the time period of Sq. Recently such nonclassical superposition of 
multiple coherent states in a Kerr medium has been experimentally achieved [45] . Formation of cat-like states 
in a finite dimensional bosonic system that admits applying a displacement operator on its ground state has 
also been studied [30] in a Kerr medium. 

In our model we study the emergence of these transitory ‘kitten’ states using the VF-distribution and its 
smoothed analog the Q-function. In the strong coupling limit A < 0.05 a;, and at specific times given by the 
rational submultiples of Tjong : {Tp,q = (p/?) TiongKp, <?) = IjP < q} density matrices comprising of a finite 
number of macroscopic coherent states with uniform angular separation on the phase space are observed (Fig. 


1^. Starting with the initial hybrid Bell state (2.9) of the composite system, the evolution of Sq in the long range 
quasi periodic regime shows (Fig. ^(b)) the existence of the local minima at rational sub multiples of Tiong- The 
presence of numerous time scales due to the qubit-oscillator interaction in the present model, however, introduces 
another novel interference related feature. In particular, the interaction-dependent linear mode with frequency 
0(a:A) causes an energy transfer, in a short time scale, between the qubit and the oscillator degrees of freedom. 
In the vicinity of the said times Tp g the oscillations of the period Tiong produce a locally minimum occupation 
on the phase space, whereas the short time period fluctuations engineer the spread of the occupation by splitting 
of the Gaussian peaks. This manifests as an ordered bifurcation {a la Figs. Ill(a2, a 4 ) and (ag, as), say) of the q 
quasi-probability peaks to 2q peaks representing mixed state oscillator density matrices, while evolving from the 
local minima to the maxima of the short time period oscillations of the Wehrl entropy. These local variations of 
Sq in the neighborhood of {p/q) Tiong are given in Figs. (ai)-(ci), where we fix p = 1; g = I, 2,3, respectively. 
Moreover, the short range time period (Tghort) associated with the splitting and subsequent rejoining of the 
peaks at a particular rational fraction (p/?) Tiong scales inversely with q. For instance, from the Figs. |4](ai)- 
(as), respectively, we observe that Tshort(Tiong) « 6750, Tshort(Tiong/2) « 3375, Tshort(Tiong/3) « 2250. This 
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suggests the scaling relation Tshort(Tiong/9) ~ (1/?) Tshort(21ong)- In other words, due to the complex nature 
of the qubit-oscillator interaction in the strong coupling regime A ~ 0.05 a; the exchange of energy is realized 
between multiple interaction-dependent modes, and effectively the high frequency quantum oscillation 0{xA) 
is frequency modulated by the low frequency component 0{x'^A). It is worth mentioning that the oscillator at 
the local minima of the Wehrl entropy (Figs. Saa, a 4 ),(b 2 , b 4 ), (ca, C4)) is close to pure states. Their respective 
von Neumann entropy S read 0.26616,0.17983, 0.53918, whereas the corresponding maxima (Figs, [^ag, a 5 ),(b 3 , 
bs), (c 3 , C5)) describe almost maximally mixed states. 

Lastly, we note the geometry of the domain on the phase space that supports the IF-distribution (Figs. 
(aa)-(ca), (a3)-(c3)). The interference pattern realized between two Gaussian peaks occurs in the intermediate 
phase space giving rise to oscillations in a direction perpendicular to the line joining the peaks. Alternate lines 
with the positive and the negative values of the VF-distribution appear with relative phase differences of tt. 
As the number of peaks increase, the interference pattern becomes more complex while being restricted within 
a regular polygon with peaks lying at its corners. The Gaussian peaks of the Q-functions (Figs. 0 (a4)-(C4), 
(a 5 )-(c 5 )) appear as smoothed versions of the IT-distributions. 



Figure 4: Transient formations of ‘kitten’ states are observed for the parameters A = 0.8 a, A = 0.01 uj,a = 2.5. 
The rows ({a}, {b}, {c}) describe the results at times Tiong, 2jong/2, Tiong/3, respectively. The column (ai, bi. 
Cl) marks the short time period oscillations (frequency 0{xA)) of the Wehrl entropy Sq. The respective times 
corresponding to the local minima in (ai, bi, Ci) are 2.9590 x 10®, 1.4778 x 10®,9.9080 x 10®, whereas the times 
of the subsequent local maxima, in turn, read 2.9612 x 10®, 1.4794 x 10®, 9.9195 x 10®. Between the minima and 
the maxima of the short time period oscillations, a doubling of the number of ‘kittens’ is observed. The columns 
(aa, ba, Ca) and (a3, b3, C3) specify the Wigner kF-distribution at the minima and the maxima of the said short 
time period oscillation, respectively. Gorresponding results for the smoothed Q-function are subsequently given 
in the columns (a4, b4, C4) and (as, bs, C5). 
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IV.C Wigner entropy and negativity 
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Figure 5: The time evolution of the Wigner entropy Sw (blue), Wehrl entropy Sq (red), and the negativity 
S\Y (black) obtained for the parameters A = 0.5w,q; = 2, where the coupling strengths read: (a) A = 0.1 w, (b) 
A = 0.2 uj. (c) The W-distribution in the stochastically stabilized domain (A = 0.2 oj) is given at the scaled time 
cut = 228, when the observed negativity remains prominent: 5w = 0.9813. 


It is also of interest to study the quantum entropy based on the modulus of the Wigner distribution \ W{j3,j3*)\ 
[29j which is a nonnegative quantity: 


Sw = -I |VF(/l,r)l log|IF(/3,r)|d'/3. 


(4.4) 


As the IF-distribution contains more information on the phase space structure of a quantum state than its 


smoothed analog the Q-function, a comparative study of the Wigner entropy Sw (4.41 and the Wehrl entropy 
Sq (4.1) is expected to throw a light on the nonclassicality of the state. It is evident from the Figs. [^(a), 
(b) that the time evolution of the Wigner entropy (4.4) and the negativity parameter (3.15) have close kinship 
with each other. This was observed for certain oscillator density functions in [5S]. In this sense the Wigner 
entropy Sw reveals the extent of nonclassicality of a quantum density matrix. We distinguish between two 
possible scenarios depending upon the qubit-oscillator coupling strength, (i) In the strong coupling regime 
(A/w < 0.1) we observe (Fig. [^(a)) a periodic structure that may be identified with the revival and collapse of 
the qubit density matrix elements reflecting the energy exchanges between the qubit and the oscillator degrees 
of freedom. The collapse, say, of the qubit density matrix elements coincides with a wider spread of the phase 
space distributions of the oscillator. The order of the time period of the revival and collapse of the qubit 
matrix elements, obtained via retaining up to the linear terms in the Laguerre polynomials, is 0(2tt/ xA). 
Therefore the variables such as the Wigner entropy Sw, the negativity Sw, and the Wehrl entropy Sq display 
similar periodic patterns in the said time scale. For a dominant value of Sw the quantum interference effects 
are overwhelming, and we, expectedly, find Sw > Sq, as an increased negativity necessitates an increment in 
the magnitude \W{/3,/3*)\ for maintaining the normalization property (3.7). This, in turn, leads to increased 
value of the entropy Sw- On the other hand, for a low negativity domain Sw 'C 1 the inequality is reversed: 
Sw < Sq. The underlying reason is that the Q-function is obtained from the W-distribution ( |3.20[ ) after 
suitable smearing with a positive definite Gaussian kernel on the phase space, and therefore it incorporates less 
information on the quantum state than the latter [17]. (ii) In the ultra-strong coupling regime (A/w ^ 0.1) all 
modes for the qubit-oscillator interaction with incommensurate frequencies are excited and a fully randomized 
interference pattern evolves very quickly. Smearing the clear periodic structures observed earlier these large 
number of interaction-generated modes lead to quasi stationary values of the phase space observables (Fig. 
(b)). However, despite the statistical stabilization of the occupation on the phase space the nonclassicality of 
the state remains prominent due to the interferences occurring between multiple modes. These interferences 
necessarily develop (Fig. i( c)) significant domains on the phase space with negative values of the W-distribution 
leading to dominant values of the negativity parameter Sw- The average value of Sw increases with that of 
the coupling strength as more interfering modes come into existence. As mentioned before, this results in an 
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(a) (b) (c) (d) (e) 



Figure 6: (a) The time evolution of the quadrature variance Vg=o for the parametric choice A = w, A = 0.1 w, a = 
0.5. The horizontal red line represents the limit Vg = 0.25. (b), (c) The plots of the VF-distribution at the 
scaled times ut — 0 and ujt = 28.80, respectively, (d) The red dotted circle indicates classical limit of the 
variance Vg = 0.25. The polar plots for the variance Vg at wt = 0 (green) and wt = 28.80 (red) are denoted. The 
least value of the variance at ujt = 28.80 equals 0.15967 and is observed at 9 = 179.55°. (e) The time-averaged 
variance (V6i=o)timeav. is plotted w.r.t. A for A = uj,a = 0.5. 


increment of the entropy Sw- In the quasi stationary state the negativity Sw is stochastically preserved. The 
stochastic stabilization of 5w occurs after a suitable decoherence time. Therefore in Fig. (b) we observe that 
except for a brief initial period the Wigner entropy Sw is consistently more than the Wehrl entropy Sq, even 
though the nonnegative Q-function may be viewed as a smeared form of the VF-distribution. Our results (Figs. 

(a), (b)) suggest that when the negativity 5w assumes more than a threshold value 5w ^ 0.4 the quantum 
fluctuations ensure the entropy relation: Sw > Sq. 

Another feature revealed in Figs. (a) and (b) is that the relative fluctuations of the Wigner entropy Sw 
overwhelms that of the Wehrl entropy Sq : I^^'q/S'qI ^ \5Sw/Sw\- It signifies that the quantum interferences 
in the evolution of the IF-distribution induce rapid reversal of its sign, whereas the fluctuations of different 
modes are significantly attenuated in the smoothing introduced towards obtaining the Q-function. 


IV.D Evolution to squeezed states 

A feature of nonclassicality such as squeezing is observed during the time evolution of the state due to the 
quantum interferences between various modes. At relatively small value of the phase space separation (a < 1) 
and for the strong coupling regime (A ~ 0.1 w), the squeezing is noticed both for the large detuning (Fig. and 
the resonant frequencies. The plot of the bF-distribution (Fig. |^(c)) makes the quadrature squeezing evident. 
The signature of the squeezing is observed (Fig. (d)) when the variance Vg of the quadrature variable, say 
at d = 0, is rendered less than its classical value 1/4. It follows from the polar plot (Fig. (d)) that at the 
scaled time wt = 28.80 the quadrature variance Vg reaches a minimum value 0.15967 at an angle 9 = 179.55°. 
The polar angle at which the minimum quadrature variance is realized varies with time and depends on the 
dynamical state of interference of the quantum modes. A qualitative understanding of the squeezing of the 
state may be described as follows. Passing to the interaction picture for the Hamiltonian (2.1), an effective 
Hamiltonian may be obtained d la |43| in an order by order perturbation theory. This effective Hamiltonian 
contains two photon terms (a^,al^) at the order 0{x) of the coupling strength. These two photon terms give 
rise to the squeezing of the state. An increase in the coupling A first enhances the squeezing as it augments the 
strength of the two photon terms in the effective Hamiltonian. Multiple photon terms, however, soon appear 
in the effective Hamiltonian [43] with increased value of the coupling strength. The resulting randomness of 
the phase relationships of the higher order terms eliminates the squeezing property in the ultra-strong coupling 
limit. It is obvious in Fig. (a) that the quadrature squeezing, even though it may be present during part of 
the dynamical evolution of the oscillator state, is not realized throughout the oscillatory cycle. The quadrature 
variance remains above the threshold value: Vg-o >1/4 during the major part of the evolution of the state. The 
instances of squeezing decreases with increasing A. To illustrate the feature we plot (Fig. e)) the time-averaged 
quadrature variance (V/=o)timeav. w.r.t. the coupling strength A. We observe that, for the parametric regime 
studied here w > A, the time-averaged quadrature variance has no squeezing property: (V/=o)timeav. > 1/4, 
and it smoothly increases with the rising A. The effective Hamiltonian approach [35], however, suggests that 
the two photon terms (a^, ^) underlying the squeezing of the state survive the smearing of the fast oscillatory 
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terms only in the limit of the large qubit frequency: w <?; A. We will return to this topic somewhere else. 

It is interesting to note that the squeezed state described in Fig. represents almost pure state of the 
oscillator at the given time as its von Neumann entropy S'|(^t= 28.80 = 0.156 90 is much less than its maximal 
value. The almost pure state of the oscillator is reciprocated by the corresponding nearly pure state of the 


qubit. For the sake of completeness we include the qubit density matix (2.14) at the given time: 


0.90760 

-0.10668-i 0.19300 


-0.10668-k i 0.19300\ 
0.09240 ) ■ 


(4.5) 


Its eigenvalues and corresponding eigenvectors read: Ai = 0.96342, A 2 = 0.03658, |Ai) = (0.46897—i 0.84843) |1)- 
0.24543 | —1) and IA 2 ) = (—0.11873 -I- i 0.21480) |1) — 0.969411 —1). The magnitude of the largest eigenvalue of 
the qubit density matix ( |2.14 1 may be regarded as measure of the purity of the state. For smaller values of 
the coherent state amplitude a < 0.1 the purity of the state may exceed 99%. Appearance of a squeezed state 
requires emergence of appropriate phase relations between various quantum modes that is facilitated by the 
proximity to a pure state. A statistical mixture of a large number of pure states is likely to destroy the phase 
relationships and, consequently, the squeezing property is eliminated. 


IV.E Mandel parameter and nonclassicality 



Figure 7: The evolution of the Mandel parameter Qm for the choice A = 0.5 w, a = 2, where the coupling 
strengths are given by (a) A = 0.01 w, (b) A = 0.08 w, (c) A = 0.15 a;. (d) The time-averaged value of the Mandel 
parameter {Qm) time w. plotted w.r.t. the coupling strength A for A = 0.5 a;, a = 2. 


In experiments allowing a direct detection of photons the Mandel parameter [30] 


Qm = {{Anf)/{n) - 1 


(4.6) 


is a convenient tool for studying the classical-quantum boundary. For the coherent state the number operator 
follows the Poissonian statistics with its signature property Qm = 0 that is regarded as the threshold for the 
classical characteristics. When the quantum correlations in the system suppress the fluctuations in the photon 
number, it assumes negative values: Qm < 0 and captures the sub-Poissonian behavior of the photon statistics 
that acts as a measure of nonclassicality. The mean value (3.26) and the variance (3.27) of the photon number 
operator obtained in our model now provide a direct construction of the Mandel parameter Qm- For the coupling 
strength A < 0.1 lv transient sub-Poissonian photon statistics is manifest (Fig. during parts of the evolution of 
the system. The time evolution of the Mandel parameter Qm in Fig. [^(a), however, passes alternately between 
the classical and the nonclassical regimes. The said sub-Poissonian behavior within short sample intervals is 
offset by the super-Poissonian photon statistics Qm > 0 during most of the time evolution of the oscillator so 
that the averaged value over a time scale tavg ^ 27r/(a;A) characterising the qubit-oscillator energy fluctuations 
always remains nonnegative. The classicality of the photon statistics, therefore, emerge via the time-averaging 
procedure even though at a shorter time scale the underlying quantum nature of the photon emission process 
is evident. At higher coupling strengths Fig. (b), (c) the randomization of the phases generated by a large 


number of incommensurate modes with the frequencies {0(a:"'A)|n = 1,2,...} partially erase the quantum 
correlations, while making the negativity of Qm less common. The effective Hamiltonian approach |43| . as 
noted in Subsec. [iwg produce multiple photon terms that give rise to the negativity of Qm- Such terms, 
however, contribute after the smoothing of the rapidly fluctuating components only in the limit A ^ uj. The 
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time-averaged properties the Mandel parameter Qm is summarized in Fig. [^(d), where, in the coupling regime 
A < 0.1 oj the time-averaged value (QM)time av. ~ 0 is maintained so that the overall photon statistics effectively 
remains to be Poissonian. After a brief transition zone the increasing coupling strength (A 3> 0.1 w) triggers 
a full randomization of the phase relations among the interfering modes, which now constitute an effective 
bath leading to progressive emergence of classical properties. The resultant stochasticity introduces a scaling 
behavior for the averaged value of the Mandel parameter: (Qm) time av. oc (A — Aq) with a suitable Aq that can 
be determined from the Fig. [Zl(d). 


V Conclusion 

Employing the generalized rotating wave approximation we have studied an interacting qubit-oscillator bipartite 
system for both the strong and the ultra-strong coupling domains for the choice of an initial hybrid Bell state. 
The evolution of the reduced density matrices of the qubit and the oscillator are obtained via the partial 
tracing of the complementary part of the respective degrees of freedom. On the oscillator phase space its 
density matrix furnishes the diagonal P-representation that is highly singular due to the presence of rapidly 
oscillating derivatives of the 5-function. Two successive smoothing performed by the Gaussian kernels on the 
singular P-representation produce first the Wigner kF-distribution, and then the Husimi Q-function. The quasi 
probability VF-distribution admits negative values due to the quantum interferences. Its negativity measure 
Sw marks the departure of the state from a positive definite distribution on the phase space. The nonnegative 
Q-function provides the Wehrl entropy Sq that acts as measure of delocalization on the oscillator phase space. 
Even though the kF-distribution may be thought of encapsulating more informations than the Q-function, the 
Wigner entropy Sw defined via the magnitude |kF| overwhelms the Wehrl entropy Sq whenever the negativity 
measure 5w is dominant. The presence of multiple time scales induced by the interaction introduces a novel 
feature in the generation of ‘kitten’ states. In the coupling range A < 0.05 a; the long-term quasi periodicity is 
observed to follow via the terms 0{x^) in the interaction Hamiltonian. In this domain the Q'function evolves 
at rational fractions of Tiong to a collection of uniformly separated Gaussian peaks representing the ‘kitten’ 
states. A shorter time scale 0((a:A)“^) now causes further bifurcation of the Gaussian peaks coinciding with 
the collapse of the qubit density matrix elements. In the chaotic A ^ a; regime all interaction-dependent 
modes of frequencies {0(a:"'A)|n = 0,1,...} with a randomization of their phases. The decoherence time 
may be estimated as the transient production time of Sq as it approaches its stochastic stabilization. Using 
the asymptotic behavior of the associated Laguerre polynomials {x) the decoherence time is estimated as 
proportional to y/X. Nonclassical features such as squeezing and negativity of the Mandel parameter arise due to 
appearance of multiple photon terms induced by the interaction. In the parametric regime studied here oj > A, 
such effects do not survive a suitable coarse graining process that smooths the rapidly oscillatory components 
in the fluctuations. 
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